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Abstract We describe recent attempts to extract the shear viscosity of the dilute 
Fermi gas at unitarity from experiments involving scaling flows. A scaling flow is 
a solution of the hydrodynamic equations that preserves the shape of the density 
distribution. The scaling flows that have been explored in the laboratory are the 
transverse expansion from a deformed trap ("elliptic flow"), the expansion from a 
rotating trap, and collective oscillations. We discuss advantages and disadvantages 
of the different experiments, and point to improvements of the theoretical analysis 
that are needed in order to achieve definitive results. A conservative bound based on 
the current data is that the minimum of the shear viscosity to entropy density ration 
is 77/5 <0.5h/k B . 



1 Introduction 



A cold, dilute Fermi gas of non-relativistic spin 1/2 particles interacting via a short 
range interaction tuned to infinite scattering length, commonly referred to as the 
unitary Fermi gas, provides a new paradigm for many strongly correlated quantum 
systems |[T] [2l ■ In this contribution we focus on non-equilibrium aspects of the uni- 
tary Fermi gas, in particular its shear viscosity |3). The shear viscosity of a liquid 
composed of weakly coupled quasi-particles can be estimated as 

f] = ^npl m jp, (1) 

where n is the density, p is the average momentum of the particles, and l m f p is the 
mean free path. The mean free path can be written as l m f p = \/[na) where a is 
the transport cross section. Equ. (HJ implies that the shear viscosity decreases as the 
strength of the interaction increases. In the unitary gas the cross section saturates 
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the s-wave unitarity bound a = \n /k 2 , where k is the scattering momentum, and we 
expect the shear viscosity to be unusually small. 

Danielewicz and Gyulassy pointed out that the Heisenberg uncertainty relation 
imposes a bound on the product of the average momentum and the mean free path, 
phnfp > K an d concluded that rj/n > h 0). This is not a precise statement: The 
kinetic estimate in equ. (Q3 is not valid if the mean free path is on the order of the 
mean momentum. A more precise bound has recently emerged from holographic 
dualities in string theory. In this context the natural quantity to consider is not the 
ratio rj/n, but r\/s, where s is the entropy density. Policastro, Son and Starinets 
showed that in j¥ = 4 supersymmetric QCD the strong coupling limit of rj/s is 
equal to h/ (4nks) G) . It was later shown that the strong coupling limit is universal 
in a large class of field theories, and it was conjectured that r\ js > h/(4nkB) is a 
general lower bound, valid for all fluids J6) . 

Are there any fluids in nature that attain or possibly violate the proposed bound? 
A fluid that saturates the bound has to be a quantum fluid (because 77 is on the 
order of Tis), and it has to be strongly interacting (because in a weakly interacting 
system the mean free path is large). It is also known that many of the model field 
theories that attain the bound in the strong coupling limit are scale invariant. All of 
these properties point to the unitary Fermi gas as a plausible candidate for a "perfect 
fluid". 

Almost ideal hydrodynamic flow in the unitary Fermi gas was first observed in 
0. Since then, a number of experiments have been performed that provide con- 
straints on the shear viscosity of the unitary gas (H 111 [Tol QT| E] \l3\ QU Q3) . In 
this work we will provide an overview of the hydrodynamic analysis of these exper- 
iments, and compare some of the estimates that have been obtained. We emphasize 
the uncertainties of these results, and point to improvements that need to be imple- 
mented. 



2 Scaling Flows 

We begin by studying the ideal (Eulerian) fluid dynamics of a non-relativistic gas 
in the normal phase. We will introduce dissipative effects in Sects. I3.1H3.3I In this 
contribution we will not discuss superfluid hydrodynamics. We will briefly comment 
on dissipative effects in the superfluid phase in Sect. 13. II The equations of continuity 
and of momentum conservation are given by 




(2) 




dv 



(3) 



where n is the number density, m is the mass of the atoms, v is the fluid velocity, P 
is the pressure and V is the external potential. In the unitarity limit the equation of 
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state at zero temperature is of the form 



n 5 / 3 / mT 



P (nJ) = —f[-^), (4) 



where f(y) is a universal function. We note that y = const - (T/Tp° m ), where 
jhom _ (37j 2 n ) 2 / 3 /(2m) is the Fermi temperature of a homogeneous Fermi gas. 
In the high temperature limit, y ^> 1, we have f(y) ~ y and in the low temperature 
limit f(y) ~ (37r 2 ) 2 / 3 £ /5, where the parameter ^ = 0.40(2) has been determined in 
quantum Monte Carlo calculations (T6|. Monte Carlo methods have also been used 
to determine f(y) for all values of y lfT71[T8t . The critical temperature for superflu- 
idity is T c /Tp" m ~ 0.15, corresponding to y c ~ 0.72. An alternative representation 
of the pressure is 

/W)=M 5/2 m^(£), (5) 

where g(z) is a universal function, related to f(y) by thermodynamic identities. In 
the high temperature limit g(z) ~ 2z 5 l 2 e l l z / (27f) 3 / 2 and in the low temperature limit 
g{z) ^ 2 5 / 2 /(157T 2 ^ 3 / 2 ). The density is 

nQi,T) = n^m^h f^j , h(z) = 5 -g(z) - zg' \z) . (6) 

The high and low temperature limits of the function h(z) are h(z) ~ 2z 3 / 2 e'/ : / (27r) 3 / 2 
(z > 1) and h(z) ~ 2 3 / 2 /(3n; 2 ^ 3 / 2 ) (z < 1). The equilibrium distribution n of a 
trapped atomic gas follows from the hydrostatic equation VPo = —noVV. The trap- 
ping potential is approximately harmonic 

V(x) = ^a>fxf. (7) 

Using the Gibbs-Duhem relation dP = ndjJ. + sdT together with the fact that the 
equilibrium configuration is isothermal we can write the equation of hydrostatic 
equilibrium as V/i = —VV. This implies that the equilibrium density is «o(x) = 
n(n(x),T) with 

MW = Mo-VW = /i fl-j:^) ) = (8) 

V i R i J m(0 r 

A scaling flow is a solution of the hydrodynamic equations in which the shape of 
the density distribution is preserved. Consider the ansatz n(x,t) = n(n(x,t),T(t)) 
where 

. . / x 2 y 2 z 2 xy \ 

MM = MoW(l- W - W - W -^yJ, (9) 

and T(t)/T(0) = /io(?)/Mo(0). Without loss of generality we have restricted the 
ansatz to rotations in the xy-plane. We note that the fluid remains isothermal during 
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the expansion. Scale invariance implies that properties of the fluid only depend on 
the dimensionless ratio T/ji. For any given fluid element this ratio does not change 
during the expansion. In particular, if the fluid element was in the superfluid or 
normal phase initially, it will stay in that phase throughout the expansion. 

The velocity field created by the scaling expansion in equ. (O is linear in the 
coordinates. We can write 

v(x,t) = iv (a x (t)x 2 + <Xy(t)y 2 + a z (t)z 2 + 2a(t)xy) + Q(t)zx x. (10) 

The parameters a, , a and Q. are related to the parameters R{ , R„ and Hq by the conti- 
nuity equation. Remarkably, the continuity equation is independent of the universal 
function h(z) in equ. ©. Introducing the dimensionless scale parameters 

the continuity equation can be written as 



2 
3 

2(a-Q) 2(a + Q) 

li + mi 



p, + -n{a x + a y + a z ) =0, (12) 



+a(a x + Oy)=0, (13) 



b x -b x a x -^(a + Q) = 0, (14) 
blX 2 a 

b y -b y ay-^—{a-n) = 0, (15) 

b z ~b z a- = 0, (16) 

where A = R y (0) /R x (0) = (O x / C0y - These equations can be solved directly in the case 
that there is no rotation, a(t) = 0. Then a = £2 = and 

x (b x b y b-\ _ 1 

(a ^' az)= Uvd' A= (w^' (17) 

The velocity field is a simple "Hubble flow", v= ((X x x, CCyy, O^z). Finally, we note 
that the entropy density is given by s = (m[i) 3 / 2 g'(T / ji). Since the entropy density 
has the same functional form as the particle density we conclude that, in the case of 
scaling flows, the continuity equation implies entropy conservation, 

^ + V-(vs)=0. (18) 
at 
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3 Elliptic Flow 



The simplest scaling flow is the expansion of the cloud after the trapping potential is 
removed lfT9l . Since the cloud remains isothermal the Euler equation can be derived 
using the Gibbs-Duhem relation dP = ndjj.. This implies that the equation of motion 
is independent of the universal function f(y) defined in equ. We get 

(O 2 1 

b ' = {b x byb-) 2 l 3 bi ' (19) 

The total energy of the expanding system is given by the sum of internal energy and 
kinetic energy, 



E = E in ,+E ki „ = Jd 3 x (^(x) + ^mnv 2 ^j 



(20) 




For the Fermi gas at unitarity the energy density $ is related to the pressure by 
S = \P. We find 



(21) 



where (0) is the internal energy at t = 0. Conservation of energy immediately 
follows from the equation of motion, equ. (TT9b . We note that the equation of hydro- 
static equilibrium, VP = — riVV, implies the Virial theorem (£) = (V) lEUl . where 
(V) denotes the integral of the potential energy over the trap. This means that the 
total energy of the trapped gas is Eq = 2£' m/ (0), where the factor 2 is due to the 
contribution of the potential energy. 

We are interested in an axially symmetric trap with coy = CO- = CO± and co x = Xco±, 
In this case we end up with two coupled equations for b± and b x . If A ^> 1 the 
evolution in the transverse direction is much faster and the equation for b± can be 
approximately decoupled, 

'b\ = %- (22) 

This equation has to be integrated numerically. The behavior at early and late times 
can be found analytically. We get 



1 + \ co\t 2 + 0{t A ) a x t<l, 
^+c + O(f- 1 / 3 ) co ± f»l, 



where 7=2/3 and co is a constant that can be determined by matching the early and 
late time behavior. Numerically, we find Co — — 1.3. For the longitudinal expansion 
the early time behavior is b x (t) ~ 1 + (X(Q±t)~ /2, and at late times b x (t) ~ const- 
X 2 co ± t. 
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The signature effect of hydrodynamics is that transverse pressure gradients cause 
the transverse radius to expand much faster than the longitudinal radius. This means 
that the two radii will eventually cross. This happens at a time 

t ao » = — (l + 0{k)). (24) 
(Ox 

We note that the crossing time only depends on the trap parameters, and is indepen- 
dent of the initial energy or the number of particles. We also note that at t ~ t cross 
the expansion is still two-dimensional, that means the volume of the system grows 
as vol ~ f 2 . The expansion becomes three-dimensional, vol ~ f 3 , at tjj ~ (A 2 fi)j_) -1 . 



3.1 Energy dissipation 

We wish to understand how the expansion is affected by dissipation. The energy 
momentum tensor of a dissipative fluid is Ily = Pdjj + mnvjVj + <5n, ; - with 

81k] = 1 faivj + VjVi - 2 - 5 l7 V • + (Sij ( V • y) . (25) 

The energy current is jf = vi(w+ jtnnv 2 ) + 8 jf with w = S +P and 5 jf = dlljjVj ~ 
KVjT. The unitary gas is scale invariant and £ = EH . Also, for an isentropic 
scaling expansion the temperature remains independent of position, and there is no 
contribution from the thermal conductivity K. We will therefore concentrate on the 
role of shear viscosity. 

Since the shear viscosity is small, we can take it into account perturbatively. The 
simplest idea it compute the amount of kinetic energy that is converted to heat. We 
have 



E = -^Jd 3 xri (v iVj + - | Sij V • v 



(26) 

For the scaling expansion given in equ. (TTOb the result is particularly simple. We get 



3U b x 



2 



d 3 xri{x). (27) 



The total energy dissipated is given by the integral of equ. d27l i over time. We first 
show that the spatial integral over r) (x) does not depend on time. In the local density 
approximation r) (x) = rj (jj.(x) ,T) . Scale invariance implies that 

r 1 (li,T)=n{n,T)a n (^j , (28) 



where a„ (z) is a universal function, and we have set Ti = 1 . In order to compare with 
the string theory bound it is also useful to define rj(ji,T) = s(fl,T)a s {T / fi), where 
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we have also set kg — 1 . We can write 

Jd 3 xri(x)^N(a n ) 1 (29) 

where 

(oc„) = -J- / d 3 xn(x,t)a„ ( T } t \ ) = t; f d 3 xriQ(x)a n ( J° . ) (30) 



is an average of a„ over the initial density distribution. Analogously, we can write 
the integral over rj(x) as S(a s ), where S is the total entropy and (<x s ) is an average 
of a s over the initial entropy density. 

The time integral over (bj_/bj_ — b x /b x ) 2 is dominated by the regime co±t ~ 1 
and converges rapidly - the integral reaches 80% of its asymptotic value at t^ ss ~ 
5 .9 ©J 1 . In the limit A < 1 we can neglect the contribution from b x . On dimensional 
grounds the integral over (bj_/bj_) 2 must be proportional to ©X- The constant of 
proportionality can be determined numerically. We find 



i: 



L ^ 2 



dt [ y- ] =0.87 o ± . (31) 



We can now compute the ratio AE/Ej nt of the dissipated energy to the initial internal 
energy of the system. In order to express the result in terms of experimentally mea- 
sured quantities it is useful to introduce the energy Ep = N£f where £p = a>(3N) l ' 3 
is the Fermi energy of the trapped gas and o = (OvO, ©-) 1 / 3 . We find 

AE 8 

--■ 0.87 J3= -2.32 j3 (32) 



E int (0) 

where the parameter j3 is defined given by 

(««) 1 («v> (S/N) 



(3Att)V3 {Eq/Ef) (WXy/i (Eq/E f ) 



(33) 



Dissipation slows down the transverse expansion of the system. For ((Q±t) 2> 1 
we have (8b±/b±) = (AE/E)/2 and, up to terms that are higher order in A, the 
change in the crossing time is directly related to the change in the expansion rate, 
(dt/t) cross = (5ix/Aj.). 

The thermodynamic quantities S/N and Eq/Ef as a function of T /Tf were de- 
termined experimentally in 11221 . Just above the critical temperature S/N ~ 2.2 and 
Eq/Ef — 0.83. The double ratio [(S/N) / (Eq/Ef)] is only weakly dependent on T, 
changing by less than 15% between T c and 4T C . In the flow experiment carried out 
by O'Hara et al. Q the cloud contained = 2 • 10 5 atoms and the asymmetry pa- 
rameter was X = 0.045. The predicted sensitivity of the crossing time to dissipative 
effects is 




(34) 

For (a s ) = 1/ (47r) this is at the limit of what can be resolved experimentally, but 
for (a s ) = 0.5 the effect reaches about 5%. An example is shown in Fig. Q] The 
solid lines show the solution of the Euler equation $1% , and the dashed lines show 
a solution of the Navier-Stokes equation (see Sect. 13.21 with (cc s ) = 0.5. The main 
effect of shear viscosity is a suppression of the transverse expansion of the system. 
We find (8t/t) cross = 6.5%, in fairly good agreement with the estimate (8t/t) cross = 
5% from equ. (f34t . 

The best fit to the data is provided by ideal hydrodynamics with (a s ) = 0. This 
is probably related to the fact that the data were taken significantly below T c , at 
T /Tf = 0.13 ± 0.05. In this regime the system is described by two-fluid hydrody- 
namics. The superfluid component has no shear viscosity but the viscosity of the 
normal component becomes very large as T /Tf — > 1231 . In a finite system, how- 
ever, the large viscosity of the normal phase is likely to be suppressed by relaxation 
time effects, see Sect. 13.51 As a consequence one observes perfect superfluid hy- 
drodynamics. The data in Fig.Q]show some deviations from hydrodynamics at very 
early and very late times. Discrepancies at early times are probably related to exper- 
imental resolution 0, while the differences at late times may be connected to the 
breakdown of hydrodynamics in the late stages of the expansion. 

We can also compute the amount of entropy generated by dissipative effects. 
Using dS = dQ/T we find 
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(On) 1 j ^ 

N 3 (3NXyP (7b/7» S K ' 

with 

I s =<OJ 1 J T dtbl 2/3 (b ± ) 2 . (36) 

For T ~ tdiss we find I s ~ 2.6 and the produced entropy is small, (AS/N) ~ 0.27 
for the conditions given above. However, the integral diverges as 7 S ~ (coxt) 1 / 3 for 
T — > °°. This result is not reliable since we expect hydrodynamics to break down at 
late times, see Sect. 13.41 



3.2 Moments of the Navier-Stokes equation 

It is clearly desirable to study the role of dissipation more directly by solving the 
Navier-Stokes equation. The Navier-Stokes equation differs from the Euler equation 
by an extra term on the right hand side, 

mn + (v • V) v ( J = -V,P - V; Sllij. (37) 

We will assume that the viscosity is small, so that derivatives with respect to ther- 
modynamic variables can be computed at constant entropy. We will also assume 
that the entropy conservation equation, equ. ( fT8l . is not modified. Physically, this 
implies that we assume that there is a reservoir that removes the heat generated by 
dissipative effects. In this case, the only correction to the equations of hydrodynam- 
ics is the viscous force in the Navier-Stokes equation. 

In general the inclusion of the Navier-Stokes term will break the simple scaling 
form of the flow. The Navier-Stokes equation also depends on the functional form 
of the pressure and the viscosity, that means we have to specify the functions f(y) in 
equ. (|4|i and a n (z) in equ. ( f28b . A simple approach that avoids extensive numerical 
work as well as model assumptions about f(y) and a„ (z) is to take moments of the 
Navier-Stokes equation. Consider the linear moments 

m J d 3 xx k n(x) (^i + ^-f^ v^j = - J d 3 xx k (VjP + Vy c5n ;; ) , (38) 

with k = 1,2,3. Since the velocity field is linear in the coordinates we find that 
the ideal fluid terms involve second moments of the density. These moments are 
related to the potential energy in a harmonic trap and, by the virial theorem, to the 
total energy of the system. The Navier-Stokes term can be integrated by parts and 
is proportional to the integral over r)(x). As a consequence, the first moment of the 
Navier-Stokes equation depends only on the parameter j3 defined in equ. d33l . We 
get 
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0.05 0.10 0.15 0.20 
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Fig. 2 This figure show the ratio (AE)/Ej n , of the dissipated energy to the initial internal energy 
as a function of the parameter /? defined in equ. B3t . The dots show the result of a numerical 
solution of the Navier-Stokes equation <39l|40t in the limit (O z /co ± — > 0, and the line shows the 
estimate given in equ. i32\ . 



(o\ 2/3 co ± ( b± b x 



{b 2 ,b x ) 2 /3b ± b L \b L b x 



(39) 



bx ~ {b\b x )^bj b x W bj- (4U) 

These equations of motion are consistent with the result in the previous section. We 
can compute the amount of energy dissipated from equ. d2Tb and ( 1391401 1. We find 

E = -\^)(^- h A 2 ■ (41) 



3 r ,mK J \b x b x 

We note that b±(t ) and b z (t) are solutions of the Navier-Stokes equation and have an 
implicit dependence on /3. As long as this dependence is smooth, bi(t,/5) bj(t,0) 
as j3 —> 0, equ. (|4TT > reduces to equ. ( T27l i at leading order in j3 . Since typical values 
of [} are quite small, we expect the estimates in the previous section to be very 
accurate. This is studied in more detail in Fig. [2] We observe that the dissipated 
energy (AE)/E is very linear in j3 even for values of (AE)/E as large as 25%. We 
note that because of turbulence solutions of the Navier-Stokes equation do not in 
general approach solutions of the Euler equation in the limit that the shear viscosity 
goes to zero. Turbulence is not present in our analysis because we do not consider 
small fluctuations. We also note that there is no continuous forcing in the case of an 
expanding gas and it is not clear whether turbulence can develop even if fluctuations 
are included. We will estimate the Reynolds number of the flow in Sect. 
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Fig. 3 Ratio 1} jn as a function of y = (inT) / « 2 / 3 for the model defined in equ. <42l The two curves 
correspond to (from bottom to top) r/j = l/(4;r),2/(4jr) with r/ = \5/(T>2\fn). The dashed line 
shows the contribution from Tji, which is the term that contributes directly to the Navier-Stokes 
equation, and the dotted line is the contribution from r/o. Note that the critical point for the onset 
of superfluidity is y c ~ 0.72. 




3.3 Scaling solution of the Navier-Stokes equation 

In this section we discuss a specific model for the density dependence of the shear 
viscosity that preserves the scaling nature of the flow even if the viscosity is not 
zero. This model allows to compute the local amount of heat that is generated by 
dissipation, and to understand some of the shortcomings of the method discussed in 
Sections [3 . 1 1 and l3~2l Consider 

P(n,T) 

77(«,r) = 77o(mr) 3 / 2 + 7 7l ^^, (42) 

where 770.1 are constants and P(n, T) is the pressure. The first term dominates in the 
low density, high temperature limit. This is the regime in which a kinetic description 
in terms of weakly coupled atoms is applicable. Kinetic theory gives l24l|25ll 

*=3^r (43) 

The second term dominates in the high density, low temperature regime. The func- 
tional form of this term is not motivated by kinetic theory. We note, however, that 
77/71 has a minimum as a function of T, as expected on theoretical l23l and phe- 
nomenological grounds |3). 

The model given in equ. d42l has two remarkable features: First, the 770 term 
does not contribute to the Navier-Stokes equation at all. The Navier-Stokes term 
V 7 [77o(mr) 3 / 2 (V 1 Vj + ■••)] vanishes since both T and V,v, are constant. Second, the 



12 



T. Schafer and C. Chafin 



r7 1 term preserves the scaling flow. Using T, V,v; ~ const we see that Vj[r]iP(n, 7)/7 
(V,-v ; - + . . .)] scales like the contribution from the pressure of an ideal fluid, V,P(«,7). 
We get 



a>| 2r/i0)| ( b ± b x 



(b 2 ,b x ) 2 / 3 b x 3T Q b ± \b ± b x 



(44) 



x {b\b x )^b x + 3T b x U± *J' 

We observe that these equations are identical to the moment equations (139140b with 
P = ri\(Q±/(3To). This is not a surprise - the Tji contribution to 77 (n, 7) vanishes 
as n —> and the assumptions underlying the moment method are satisfied. The Tjo 
term, on the other hand, does not vanish as n — > 0, and it cannot be included in the 
moment equations (it makes an infinite contribution to the integral over 77 (x)). 
Using the identification j3 = 7]i©x/ (37b) we can write 

771 1 (46) 



3(3AA0V3 (7b/7» ' 

which shows that any bound on (a n ) obtained using the methods of Sect. I3.2l can 
be translated into an estimate of Tji, r\\ = 3(To/EQ)(a„) . Near T c this implies that 
rii ~ 0.76(a„). We note that the relation between T7 1 and (a„) is precisely what one 
obtains if the trap average of 77 (x) is computed from the rji -term only. The situation 
is more complicated if the contribution from 770 is taken into account. The ratio T] jn 
is given by 

^ = W /a +%W (47) 

n y 

with y = (mT)/n 2 ' 3 . Since /(0) = const and f(y) ~ y for y ^> 1 this function has 
a minimum, see Fig. [3] The figure also shows that (17 /«)„„„ receives significant 
contributions from Tjo- It is clearly unsatisfactory that our analysis has no sensitivity 
to this term. We will return to this issue in Sect. 13.51 

Using the explicit form of 77 (n, 7) we can also address the question where the 
energy is being dissipated and how much reheating is taking place. We first consider 
the contribution from 771. The energy dissipated is 



477! fb ± \ 2 P(n,T) 



3 \b ± 



(48) 



For a Fermi gas at unitarity the energy density is related to the pressure by S («, 7) = 
(3 /2)P(n, 7). Equ. (08]) implies that the energy dissipated is proportional to the local 
internal energy density. The source of the dissipated energy is the reduction in the 
kinetic energy density relative to its value in ideal hydrodynamics. The local kinetic 
energy density is 

^n = \n(^)A. (49) 



Scaling Flows and Dissipation in the Dilute Fermi Gas at Unitarity 



13 



Since the kinetic energy density differs from the spatial distribution of the dissipated 
energy there has to be a dissipative contribution to the energy current. This current 
is given by 8j £ = (0,8j*,Sj*) with 

8 Jz = Vz 8n a = : y/ . (- J , (50) 

and Sj £ = 8 j£(z O y). The dissipative current flows from the outer edge of the 
cloud, where the kinetic energy is peaked, to the center of the cloud, where the 
pressure is largest. 

Energy dissipation leads to reheating. The change in temperature is AT = 
(A$)/cy. The time evolution of the temperature is governed by 

t-^ikhm 1 * 

where the first term is related to the adiabatic expansion of the system, and the 
second term is the dissipative correction. Note that if cy ~ S/T, which is the case 
in the high temperature limit, then reheating will preserve the fact that the cloud 
is isothermal. In general the behavior of the specific heat is more complicated and 
dissipation produces a temperature gradient. The relative importance of reheating is 
governed by the parameter (rji(0±/To)(P/ (cyT)). In the high temperature limit we 
can use P ^cyT and this expression reduces to the parameter j3 defined in equ. (|4oT ). 
Reheating becomes important at a time co±t ~ j3~ 3 . Since p is typically very small, 
this occurs very late during the evolution of the system. 

A similar analysis of the effects of T]o leads to a number of puzzles. The energy 
dissipated is independent of density, and the total energy dissipated over all space 
is infinite. There is no change in the kinetic energy, and the source of the dissipated 
energy is the viscous correction to the energy current. This current flows into the 
system from spatial infinity. The relative importance of reheating is governed by the 
parameter (rio(0±/T())((mT) 3 / 2 /n), which is always large in the dilute region of the 
cloud. 



3.4 Breakdown of hydrodynamics 

The constant term 77 ~ rjo(mT) 3 ^ 2 in the shear viscosity dominates in the dilute 
outer regions of the cloud, and the difficulty in understanding the effects of this 
term must be related to the breakdown of hydrodynamics in the dilute regime. A 
standard criterion for the applicability of hydrodynamics is the condition that the 
Knudsen number Kn = l m f p /L, the ratio of the mean free path to the system size, is 
much less than one. In the dilute regime the mean free path is given by 
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1 3 mT 

Imfp = — = -. ■ (52) 

no An n 

The density is given by equ. ©. In the dilute regime we can use the high temperature 
limit of h(z), but the scaling arguments in the following are independent of the 
functional form of h(z). For a comoving observer the density scales as n ~ (m/z) 3 / 2 , 
and the mean free path scales as l,„f p ~ ^/(m 1 / 2 ^ 3 ' 2 ). The evolution of T and jx 
is governed by the scaling relations discussed in Sect. [2] We may use, in particular, 
that T I \l ~ const and fi ~ jx(0) / (b\b x ) 2 ' 3 . We conclude that in a comoving fluid 
cell 

7 / , \ 1/3 

hnfp I "x 



During the two-dimensional expansion the Knudsen number is dropping, which 
implies that the hydrodynamic description is becoming more accurate. In the late, 
three-dimensional stage, the Knudsen number is constant. 

A more accurate criterion can be obtained by using a characteristic length or time 
scale derived from the flow profile. Hydrodynamics is based on a derivative expan- 
sion of the energy momentum tensor, and the validity of hydrodynamics requires 
that dTIij is small compared to the ideal fluid stress tensor. Consider the ratio of 
the moments of the ideal and dissipative terms on the RHS of the Navier-Stokes 
equation 

fa V * P > = ( p ) (54) 

where (.) denotes an integral over d 3 x and the index k is fixed. The ratio (77 /P) (V • v) 
has a simple interpretation in kinetic theory. For a dilute gas 77 ~ npl m j p ~ pw 2 T m # 
and P ~ pu 2 , where n is the particle density, p is the mass density, p is the average 
quasi-particle momentum, u the average velocity, and t m ft the mean free time. The 
ratio V ■ v ~ T^ ; , defines a characteristic expansion time. The quantity 

5(V-v)~5* (55) 

' *exp 

measures the ratio of the mean free time over the expansion time. Hydrodynamics is 
valid if T m fi <C T exp . We observe that for 77 ~ P the freezeout criterion is independent 
of position and only a function of time. We get 

\ (V.VJ - I = pjL^M"' , (56) 

where we have assumed that the expansion is two-dimensional. We note that the 
relevant parameter is the quantity j3 defined in equ. d46l i. Freezeout occurs at 
(ffljfo) ~ /3~ 3 . For typical values of j3 we find that tp 3> t cmss S> t^, where 
tcmss ~ (Cx^-) 1 is the crossing time, and t c n ss ~ 5.90J 1 is the characteristic time 
for dissipative effects. 
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The freezeout time defined by equ. d56l l is very long, and the physical freezeout 
is determined by the viscous effects in the dilute part of the cloud. In the case of a 
spatially constant shear viscosity we find 

P y zv P \b ± J 8V2(3XNy/ 3 x \ibjRjj 

where we have used P =nT as well as the low density (high temperature) limit 
of ito(x), see equ. ©. The radius parameter Rj is defined as Rj = 2To/(m(of). The 
condition (77 /P)(V z v z ) determines a freezeout surface Xf r (t). This surface is initially 
atx, 3> Ri, but it moves inward as time increases and reaches the origin at a time tfr ~ 
©J 1 (3XN)(Tf /Tq) 6 . This time is also parametrically very long, but the freezeout 
time at a characteristic distance x, ~ biRi is significantly smaller. 

Finally, we wish to mention one more quantity that characterizes a viscous flow. 
The Reynolds number Re is defined as the ratio of inertial and viscous forces in the 
system. In the case of a scaling flow with 77 ~ P this ratio is independent of position 
and only a function of time. We find 

Re = ^K r b ± b x ~^. (58) 

The Reynolds number is zero initially, but it grows quickly, reaching Re ~ 
at ((0±t) = 1. For typical experimental parameters /3 _1 ~ 100, which is large but 
not large enough to cause instabilities. At later times even larger values of Re are 
reached, but at these late times the system is simply free streaming. A constant 
contribution to the viscosity does not lead to a viscous force, and does not directly 
contribute to the Reynolds number. 



3.5 Relaxation time approach 

The discussion in the previous section does not fully resolve the problems caused by 
the dilute regions of the cloud. If the shear viscosity is proportional to the pressure 
then the system freezes out at some time tfr. For values of r]/P implied by the 
data this time is much larger than the characteristic time for dissipative effects in the 
evolution of the system, and the estimates in Sect. |3.1||3.3| are internally consistent. If 
the shear viscosity is constant then there is a freezeout surface which moves inward 
as a function of time. This implies that the integral in equ. d27l > and ( f38l > should 
be restricted to the region enclosed by the freezeout surface. However, in order for 
energy to be conserved, and for viscosity to have an effect on the evolution of the 
system, we would have to include an external force on the freezeout surface. 

An approach that can describe the effects of freezeout without the need to intro- 
duce an artificial surface is second order viscous hydrodynamic |26l . The second 
order formalism takes into account terms with two derivatives of the thermody- 
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namic variables in the dissipative correction to the stress tensor and energy current. 
In general, the second order formalism contains a large number of new transport 
coefficients. A phenomenological ansatz that has proven to be useful in many dif- 
ferent applications is to treat the viscous part of the stress tensor as an independent 
hydrodynamical variable which satisfies a relaxation equation 

x R j t sriij = - dn.j + 8n*j s , (59) 

where Tr is the relaxation time and Sn^ s is the Navier-Stokes expression for the 
viscous contribution to the stress tensor, equ. d25l ). An equation of this type was 
first introduced by Maxwell and Cattaneo in the context of heat transport. More 
recently, time or frequency dependent viscosities were considered in the study of 
Bose condensed gases in ||27ll28l . In relativistic hydrodynamics relaxation equations 
for the viscous stress tensor are used in order to restore causality, see the review |29l . 

Scale invariance implies that Tr(«, T) = w(mT /n 2 ^)/T where w(y) is a universal 
function. In the dilute limit y ^ 1 the function w(y) can be calculated in kinetic 
theory which gives Xr = r\/(nT) QUI . This result corresponds to the estimate for 
%nft given in equ. ( 1551 1. The relaxation equation d59l requires an initial condition for 
the viscous stress 5Hy. If is natural to assume that <5fl, ; - = at t = 0. In the center 
of the cloud Tr is small and the viscous stress quickly relaxes to the Navier-Stokes 
result. In the dilute region Tr °o and the viscous contribution to the stress tensor 
remains zero. This implies that even a spatially constant shear viscosity leads to a 
spatially varying 5fl, 7 and a non-zero drag force. This drag force is largest near 
the freezeout surface and breaks the scaling nature of the flow. This means that a 
detailed study of the Israel-Stewart equations will require numerical solutions of 
the hydrodynamic equations. We can estimate the effect of the relaxation time by 
computing the energy dissipation. We have 

£ = ~Jd 3 x Sn^ (v iVj + Vj Vj - | S u V • v\ , (60) 

where 5JTy is determined by equ. (l59l . The simplest approximation is to set 5fl i; = 
8n?j S inside the freezeout surface and SJIy = outside. 

In order to obtain more accurate estimates we have to solve the differential equa- 
tion ( 1591 . As in Sect. 13.1 I we may compute 8nlj S from the solution of ideal hydro- 
dynamics. The relaxation time can be calculated using the high temperature result 
for the density profile. We find 

457T 1 / T V 4/3 / A 4\ 

^ = ^(m^ UJ K exp [m + w ' (61) 

which has the same functional form as the freezeout criterion in equ. (1571 . The 
viscous stress tensor 5fl; 7 is determined by integrating equ. (1591 and the dissipated 
energy can be computed from equ. (160t . By comparing AE with equ. ( l27l we can 
express the result in terms of an effective (a„). This quantity is shown in Fig. [4] We 
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Fig. 4 Trap average (a,,) = (rj) /N computed from a relaxation time equation with 77 = t)o(mr) 3 ' 2 
and Tg = r]/(nT). Contrary to the pure Navier-Stokes case T# — > the ratio {rj)/N depends on the 
number of particles and the trap geometry. Here we have chosen N = 2 ■ 10 s and A = 0.045. The 
solid shows the result for the elliptic flow field, and the dashed line corresponds to the transverse 
collective mode, see Sect. [5] 



observe that (a„) grows with temperature as (a„) ~ r 3 , much faster than one would 
expect from the relation 77 ~ r 3 / 2 . 

There are no data for elliptic flow at temperatures above T c , but we will com- 
pare the relaxation time result to collective mode data in Sect. [5] We note that at low 
temperature the effective (a„) is the same for expanding and oscillating systems, but 
that at high temperature the two systems behave differently. In the expanding system 
the hydrodynamic expansion time T exp continues to increase during the expansion, 
whereas the period of the oscillation provides a fixed hydrodynamic time scale in 
the case of the collective mode. The viscous relaxation time Tr increases with tem- 
perature. This implies that for the collective mode we eventually get t/r > x exp and 
the effective (a„) starts to decrease. In the expanding system, on the other hand, the 
relaxation time can always match the expansion time and (a n ) continues to grow 
with temperature. 



4 Expansion from a rotating trap 

The expansion from a rotating trap was studied in [ 14 1 . Rotating gases are of interest 
for a number of reasons. The quenching of the moment of inertia in a superfluid Bose 
gas was used as a signature of superfluidity 1311 . The remarkable discovery in lfl4l 
is that in a Fermi gas at unitarity the suppression of the moment of inertia is also 
observed in the normal phase. It is clearly of interest to determine to what extent 
this discovery places constraints on the shear viscosity 11321 . 
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4.1 Ideal fluid dynamics 



The Euler equations for a Bose gas with P ~ n were derived in PP . The result is 
easily generalized to a Fermi gas at unitarity lfl4l . As in the case of a non-rotating 
trap the equations are independent of the temperature and the universal function 
fly) in equ. ©. We have 

2 

a x + a 2 + a 2 -Q 2 = ^ (62) 
H 

Ov + c^ + cr-^ 2 = ^-f (63) 

a- + a} = (64) 
bi 

a + a la x + Oy) = ( 65 ) 

Q + Q{a x + Oy) = 0. (66) 

These equations have to be solved together with the continuity equations ( ll2H16b . 
In all there are ten coupled equations. In the case of a rotating trap there is no initial 
expansion, 0£, (0) = 0, but either a(0) or X2(0) (or both) are non-zero. If the initial 
flow is purely irrotational then a(0) = (O rot , where co rot is the angular velocity of the 
trap. If the flow corresponds to rigid rotation then X2(0) = (O rut . Below the critical 
temperature the flow of the superfluid component must be irrotational, but above T c 
both rotational and irrotational flows are possible. 

The equations simplify in the experimentally relevant case of strongly deformed, 
slowly rotating traps, co rot < (0 X <C (0±_ with o ± = COy ~ (O z . In this limit the motion of 
the fluid is dominated by the transverse expansion of the system. Up to corrections 
of order <9(A 2 ) or O((t0 rof /C0x) 2 ) we have 



1 + i»|f 2 + 0(f 4 ) o ± f<l 
^+c + O(r- 1 /3) fi) ± r»l. 



MOM E^W"^ mTt^l < 67 > 



as in the case of a stationary trap. The orientation of the expanding cloud is described 
by the parameter a defined in equ. <fTTT >. We find 

a{t)C " {-ffif, co ^ t » 1 ( r< ^/), (68) 

where c a is a constant. Below we will show that c a = y. At very late times, f > ~ 
1 / (A 2 a>x), we find a(t) ~ 1/f 2 . The result d68l holds irrespective of the nature of 
the initial rotational flow. The parameter a(t) can be related to the angle of the cloud 
with respect to the x-axis, 
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Fig. 5 Time evolution of the angle of the major axis of a rotating expanding cloud after release 
from the trapping potential. The data are taken from 1141 . The two data sets were obtained with 
initial energies E/Ep = 0.56 and 2.1. The solid line shows the prediction of ideal fluid dynamics, 
and the dashed lines shows the solution of the Navier-Stokes equation for j8 = 0.061. Using an 
entropy per particle S/N ~ 4.8 this value of j8 implies a shear viscosity to entropy density ratio 
{a s } = 0.60 



ak 2 blbl 

At early times, ay -^i 1, the angle is proportional to the rotational frequency of the 
trap, 6 — CO rot t. The angular motion speeds up as b y X approaches b x . The angle goes 
through 45° at 

t 45 o = ^1 (70) 

(Ox 

which is the identical to the crossing time in equ. ( f24b . At late times, and up to cor- 
rections of O((0 mt /(0±), the angle approaches 90°. The velocity field is dominated 
by the transverse expansion of the system. In the limit co ro r < ®* <C (0± the velocity 
fields a,- are identical to those in the non-rotating case. We have 

[ft c ± r«l, 
\ l/t ©j_*> 1, 

and OCx = 0(X 2 ). The rotational components of the velocity field decay quickly. If 
the initial flow is irrotational, a(0) = (Q rot , then 

a(t) ~ow(l-a)|f 2 ) (72) 
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Fig. 6 Time evolution of the parameters a, b x , b y , b z that characterize the scaling expansion out of 
a rotating trap. Note that in this case cuy and (O z are not exactly equal, and that the time scale is 
different from Fig. \5\ Here, we only show the early evolution of the system. Solid lines show the 
solution of the Euler equation, and dashed lines show the solution of the Navier-Stokes equation 
for /3 = 0.077. 



for (d)±t) < 1. For ((Q±t) > 1 the rotational component of the flow is small, 
(a/co ro t) <C 1, but the remaining flow decays slowly, a ~ f _1 for t < and 
a ~ t~ 2 for t > t^. In ideal hydrodynamics an initially irrotational flow will re- 
main irrotational, £2(t) — 0, for all t. If the initial flow corresponds to rigid rotation, 
Q. (0) = aw, then the early time behavior is given by 

Q{t)~(o rot (\-^(olt 2 ^j . (73) 

An initially rigid rotating flow induces a non-zero irrotational flow. For ((0±t) > 1 
both components of the velocity field become much smaller than (Q rot . 
The angular momentum is given by 

L z = am(n(x 2 - y 1 )) + Qm(n(x 2 + y 1 )) + {a x - CCy)m(nxy), (74) 

where n is the density and (.) is an integral over the cloud. The moment of inertia 
of a rigid rotor is 7 n? = m(x 2 +y 2 ), and the irrotational moment of inertia is Ij rr = 
m(x 2 —y 2 )- We have 



™("* 2 ) = TTT 77 T7 - ( 7 5) 



X 2 b 2 
1 - WtdbJ>y¥ a* 
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Fig. 7 Time evolution of the parameters a and 
components of the velocity field. Parameters were 
of the Euler equation, and dashed lines show the 
0.077. 



Q. which control the irrotational and rotational 
chosen as in Fig. [6] Solid lines show the solution 
solution of the Navier-Stokes equation for j8 = 



^ab 2 x b 2 y Lp 

%(ab x b y ) 2 at 



m{wcy) = - . 2 , (77) 



(78) 



where the scale is set by 



6 A 2 / 3 \E F 

In the experiment of Clancy et al. (co rot /(O x ) ~ 0.4 and Lq/N ~ \3\{Eq/Ef). For 
Eo/Ef = 1, which is in the normal phase, the angular momentum per particle is 
507i. 

At early times the trap is strongly deformed and 7 ng ~ . When the cloud be- 
comes almost spherical the irrotational moment is much smaller than the rigid mo- 
ment of inertia, /,,-,- <C I r ig- However, at times ((0±t) > 1 the angular momentum is 
mainly carried by the last term in equ. d74i i. which is related to the transverse ex- 
pansion of the system. This is true irrespective of the nature of the initial rotational 
flow. For ((Oj_t) > 1 we have Oyfn(nxy) ~ (c a / y)(co mt / CO x )Lq. Angular momentum 
conservation then fixes the constant c u in equ. d68l ), c a = y. At very late time, t > t^, 
the angular momentum is shared among all the terms in equ. d74l i. and the relative 
size of the different contributions depends on the initial conditions. In practice, of 
course, hydrodynamics is no longer applicable at t > t^. 
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Fig. 8 This figure shows different contribution to the total angular momentum of the expanding 
cloud as a function of time. The angular momentum is given in units of the quantity Lq defined 
in the text. The curves labeled irrotational, rigid, and expansion show the (x 2 — y 2 ), {x 2 +y 2 }, and 
(xy) contributions. The solid and dashed lines correspond to ideal and viscous hydrodynamics, 
respectively. The solid black line shows the (conserved) total angular momentum. 



4.2 Dissipation 



The effects of dissipation on the expansion from a rotating trap can be studied in 
close analogy with Sect. I3.U3.5I The rate of energy dissipation is 



[at 



a 2 - a T Oy- a T a-- a,,a- + 3a 2 ) Jd i xr\(x). (79) 



For a x ~ (Xy ^> a z , a this expression reduces to the energy dissipated by the trans- 
verse expansion of cloud, see equ. (f2Tb . This implies that the main effect of dissipa- 
tion is to slow the transverse expansion of the cloud, and to delay the time ?45o . This 
delay is exactly the same as the delay in the crossing time in equ. d34l i. We have 



— ) 

t ; 45 o 



1/(4*) 



1.3 - 10 5 
N 



1/3 



03\ 1/3 [S/N\ ( 2.1 
A J I, 4.8 J \E /E F 



(80) 

We can confirm this estimate by solving the Navier-Stokes equation. The Navier- 
Stokes equation can be derived using the moment method described in Sect. 13.21 
As before, an equivalent set of equations can be obtained from the viscosity model 
given in equ. (l42l . We get 1331 



a x - 


- - 


-a 2 -£2 2 = 


®, 2 1 
bl 1 
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co± 


'2 

- 3 a x 


CCy- 




-a 2 -Q 2 = 


G>y I 

b 2 1 
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61 

co± 


'2 

3«> 



(a y + a z ) + -ab~a 
(a x + a z ) + \-aX 2 b 2 a 



(81) 
(82) 
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Fig. 9 This figure shows the angular velocity of the rotating cloud as a function of time. The solid 
line shows the solution of the Euler equation, and the dashed line is the solution of the Navier- 
Stokes equation for /3 = 0.077. The thin dashed line shows the result for the angular velocity 
obtained by rescaling the solution of the Euler equation by a factor 1 + (St/t)^ ~ 1.1. The dis- 
crepancy between the Navier-Stokes prediction and the rescaled Euler result in the regime where 
® is large is due to the rotational component of the flow. We note that / = L/0 is the moment of 
inertia. 



ct 7 + a, 



d + a (a x + GCy) = co 




(a x + Oy) 



b 2 x + X 2 b 2 y 



Q + Q (a x + GCy) 



3M 2 



- (a x - Oy) 



X 2 b], 



X 2 b 2 b 2 



■ a 



(85) 



These equations are independent of the functional form of the pressure. A solution 
of the Navier-Stokes equation for the trap parameters and initial conditions in lfl4l 
is shown in Fig. [5] The experimental data were taken at E /Ep = 0.56 which is in the 
superfiuid phase, and E/Ep = 2.1 which is significantly above the phase transition. 
Similar to the low temperature data for pure transverse expansion in Fig. [T]the low 
temperature result for a rotating cloud shows no dissipative effects, and the best fit 
to the data is provided by ideal fluid dynamics. 

The data for E /Ep =2.1 clearly show a delayed expansion. We find (St/t)^ ~ 
0.063. Using (8t/t) ~ 1.16/3 from equ. (O we estimate j3 ~ 0.057. This estimate 
is quite accurate, the best fit of the Navier-Stokes solution to the data is obtained 
for J3 = 0.061. Using N = 1.3 • 10 5 , A = 0.03 03 and (S/N) ~ 4.8 J22) we obtain 
(ot s ) — 0.60. The measurements were extended to values of E /Ep between 0.56 and 
2.1 in ll32l . This work reports values of r]/s as small as (a s ) — (0.0 — 0.4). Note 
that in this regime it becomes very difficult to measure the viscosity accurately. A 
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Fig. 10 Time evolution of the amplitude of the transverse breathing mode. The black line shows 
the solution of the Euler equation and the solid green line is the solution of the Navier-Stokes 
equation for j3 = 0.05. The dashed green line is the damped cosine function given in equ. 490l >. The 
trap frequency was chosen to be C0± = 1696 Hz as in 1101 . 



value of (ofj) =0.1 affects the measured angle of the cloud by less than the with of 
the lines in Fig. [5] 

A more detailed study of viscous effects on the evolution of the system is shown 
in Figs. |6][8] We observe that viscosity slows down the evolution of the scale pa- 
rameters b y ,b z and a. More interesting is the effect on the velocity fields a and £2. 
Viscosity converts a fraction of the irrotational velocity field a into the rotational 
velocity field Q. This is also seen in the breakdown of the angular momentum, see 
Fig. [8] The rotational component of L z is not large, but it does lead to an observ- 
able effect in the angular velocity of the cloud. Fig. [9] shows that viscosity leads 
to a decrease in ©. During most of the evolution this effect is dominated by the 
delayed expansion, but for t ~ ?45 there is an extra reduction which is due to an 
increase of the effective moment of inertia / = L/0 caused by the rotational flow. 
Unfortunately, the experimental data are for ©(f) are not sufficiently accurate to 
demonstrate this effect. 



5 Collective Oscillations 



In order to study collective oscillations we consider the Euler equation (fT9] l in the 
presence of the trapping potential. The equation of motion is 
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* l= (w^T**'- (86) 

The equilibrium solution is b x = b y = b- = 1 . We now consider small oscillations 
around the equilibrium, b((t) = 1 +a,e ,a ". The linearized equation of motion gives 

co 2 ai = cof (leu + rE fl i^ ' ( 8? ) 

which was derived in l34l [35l [36l using slightly different methods. For the radial 
breathing mode a y = a z = a±, a x = we get ft) 2 = 2(1 + y)o 2 = (10/3)©j\ The 
energy dissipated can be computed from equ. (l27l i. We find 



AE 



(88) 



where AE is the energy dissipated per period, E osc is the energy of the collective 
mode, and j3 is the parameter defined in equ. d33l ). We note that the amount of energy 
dissipated in one period of the transverse breathing mode is about three times larger 
than the energy dissipated by transverse expansion, see equ. (l32l . 

We can also derive a Navier-Stokes equation, either by taking moments as in 
Sect. 13.21 or by using a simple scaling form of the shear viscosity as in Sect. 13.31 
For the transverse breathing mode we find 

*x = 4-«4»x-^. (89) 
b[ b ± 

If j3 is small then this equation is approximately solved by a damped oscillating 
function. We have 

b±(t) = 1 +aj_cos(fiiM)exp(— rt). (90) 

Comparison with equ. (l88l gives F = j5co±. The main feature of collective modes is 
that the viscous term exponentiates so that even very small values of j3 are experi- 
mentally accessible. In Fig. ( TTOb we show a comparison between an exact solution of 
equ. d89t for j3 = 0.05, a±(0) = 0.25 and the approximate solution ( 190) . We observe 
that the approximate solution is extremely accurate. 

The experimentally measured damping rate can be used to estimate (<X S ). We 
have 

<*>-^ (3 (!)(!)■ 

In Fig. [TT] we show an analysis of the data obtained by Kinast et al. J9) using 
equ. j9TT l. This plot is very similar to our earlier analysis f37l (see also 11381 [39l ), 
except that the temperature calibration and thermodynamic data have been updated 
using the recent analysis published in 11221 . 

There are a number of important checks on the interpretation of the damping date 
in terms of viscous hydrodynamics that should be, or have already been, performed. 
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Fig. 11 Trap average (a s ) = (rj/s) extracted from the damping of the radial breathing mode. The 
data points were obtained using equ. (9l\ to analyze the data published by Kinast et al. (9). The 
thermodynamic quantities (S/N) and E^jEp were taken from 1221 , The solid red and blue lines 
show the expected low and high temperature limits. Both theory curves include relaxation time 
effects. The blue dashed curve is a phenomenological two-component model explained in the text. 



Viscous hydrodynamics predicts that the monopole mode in a spherical trap is not 
damped at all. This prediction is quite striking, but it has never been tested. Viscous 
hydrodynamics also predicts simple relationships between the damping constant of 
the radial breathing mode and the radial quadrupole as well as the scissors mode 
1571 . These predictions agree qualitatively with the data obtained by the Innsbruck 
group, but there are some structures in the data that do not fit a simple hydrody- 
namic description. Finally, hydrodynamics predicts that the damping rate decreases 
as N~ 1 ^. This prediction does not agree with the data published in (9). We note, 
however, that Kinast et al. only checked the scaling behavior at very low tempera- 
ture, and that relaxation time effects may modify the particle number scaling. 

We can also compare the results in Fig.[TT]to theoretical prediction for the shear 
viscosity in the low and high temperature limit. In the high temperature limit the 
viscosity is independent of density and the main source of dissipation is the finite 
relaxation time, see Sec. 13.51 In the case of periodic motion the relaxation time 
equation ( l59l is easily solved. The dissipated energy is given by equ. ( 133188b with 

(On) = 77o('«n 3/2 / d 3 x - . 1 (92) 
J 1 + (O z TR{n{x)) z 

We will use the kinetic theory result Tr(«) = 7]/{nT) with r\ = rjo(»ir) 3 / 2 . In the 
high temperature (low density) limit we can use the classical expression for the 
density profile n(x). In this case the integral over x can be done analytically. We find 
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Fig. 12 Damping rate of the radial breathing mode in units of the transverse trapping frequency. 
This figure only shows the contribution from the dilute corona, computed using the relaxation time 
approach. The solid line corresponds to N = Nq = 2 • 10 s , A = 0.045 as in (9)- The long dashed 
and short dashed lines corresponds to N = 5No and N = 0.27Yo, respectively. 



. 457T fT\\. 

where const = 1 125 • 3 1 ' 3 7T 2 /64 ~ 250.1, and Li a (x) is the polylogarithm function. 
In the limit T <C Tf the result scales as (a„) ~ y 3 log(y) 3//2 with y = T/Tp. For 
r > 7f we get (a n ) ~ y . These results imply that both the temperature scaling 
and the particle number scaling differ from naive expectations. The shear viscosity 
scales as 77 - T 3 / 2 , but (a n ) ~ T 3 log(r) 3 / 2 at low T, and (a„) ~ at high T. 
Also, the scaling of the damping rate with N is T ^ N~ 1 ^ log(Af) 3 / 2 at low T and 
r ~ A^ 1 / 3 at high T, see Fig. [12] This implies that there are temperature regions in 
which the dependence of the damping rate on is small. 

The prediction of equ. d93l is shown as the solid blue line in Fig.Qj] We observe 
that the relaxation time model agrees well with the data for T ~ (0.5 — 0.8)7>. For 
temperature less than 0.57> the observed damping rate is bigger than the prediction 
of the relaxation model. At very low temperature the shear viscosity is expected to 
be dominated by the phonon contribution 11231 

77=0.018« — . (94) 
\ ml I 

At low temperature we can compute the trap average by using the zero temperature 
profile. We find (a„) ~ 1 .5 • 10~ 5 (7> /T) 5 . This result becomes large for T/T F < 0.1. 
In this regime relaxation time effects are important, and (a„) at finite frequency goes 
to zero as T — > 0. 



const 



(AAO 2 / 3 \T F 



(93) 
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Fig. 13 Time scales relevant to the expansion of a unitary Fermi gas from a deformed trap. The 
inverse trap frequency is coj 1 = 0.024ms. The scale t acc is the characteristic time for hydrody- 
namic acceleration, where we have defined t = t acc to be the time when 80% of the initial internal 
energy has been converted to kinetic energy. The characteristic time for viscous effects, fa, is de- 
termined by the condition that the dissipated energy AE has reached 80% of its asymptotic value. 
The freezeout time tp- is quite uncertain. Here, we show the time at which, for Tg/Tp = 0.21, the 
freezeout surface reaches the point x ± = b ± R ± . The crossing time t cr is the time at which the sys- 
tem becomes spherical. The time t^d at which the expansion becomes three-dimensional is bigger 
by another factor X ~ 1 . 



Neither the low temperature nor the high temperature result provide a good 
description of the data in the regime T ~ (0.15 — 0.40)7>. The dashed blue line 
in Fig. QT| shows a purely phenomenological fit based on the functional form 
77 = 77o(mr) 3 / 2 + riin 5 / 3 /(mT) with T70 = 15/(320r) ~ 0.264 and 771 ~ 0.06. In 
this case the minimum value of r\/n is 0.24 which occurs below the phase transition 
atmr/n 2 / 3 ~0.47. 



6 Summary and Outlook 

A special feature of the hydrodynamics of a unitary Fermi gas is the existence of 
simple scaling solutions of the equations of ideal fluid dynamics. These solutions 
are independent of the equation of state, the initial temperature and the number of 
particles. The only time scales in the problem are the trap frequencies, see Fig. [13] 
The existence of scaling solutions is related to the constraints imposed by scale 
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invariance on the equation of state, and to the harmonic character of the confinement 
potential. 

The properties mentioned above make scaling flows an ideal class of solution to 
study the effects of shear viscosity. In this contribution we focused on three classes 
of experiments, expansion from a deformed trap ("elliptic flow"), expansion from 
a rotating trap, and damping of collective oscillations. These experiments provide 
somewhat complementary information, and they have different advantages and dis- 
advantages: 

• In the case of collective modes the effect of shear viscosity exponentiates, and as 
a consequence the damping of collective modes is sensitive to very small values 
of the shear viscosity. Collective modes also have the advantage that qualita- 
tively the effect of dissipation is very simple: The kinetic energy of the collective 
mode is converted to heat, so that at the end of the evolution the system is again 
stationary, but the temperature is increased. In the case of flow experiments the 
situation is more complicated. Dissipation converts kinetic energy into heat but 
unless the system freezes out first, the internal energy is eventually converted 
back to kinetic energy. Because of the second law of thermodynamics, the final 
state of viscous hydrodynamics must differ from that of ideal hydrodynamics, but 
the differences can be subtle, manifesting themselves in violations of the simple 
scaling formulas for the density and the velocity field. 

• The transverse expansion experiments provide detailed information about the 
time dependence of the density and flow profiles. This information can be used 
to understand the breakdown of hydrodynamics, for example by studying devi- 
ations from the simple linear velocity profile predicted by ideal fluid dynamics. 
Transverse flow experiments may also show a different, and possibly smaller, 
sensitivity to relaxation effects. Fig.|4]shows that, for T/Tp < 0.4, the relaxation 
time estimate of the trap averaged dissipation due to the spatially constant part of 
the shear viscosity is similar for transverse flow and transverse collective modes. 
However, the local response of a rapidly expanding cloud is likely to be different 
from that of an oscillating system. 

• The expansion of a rotating cloud is sensitive to a new viscous effect, the con- 
version of an irrotational flow v ~ V(xy) to a rotational flow v~2xx Contrary 
to the slowdown of the transverse expansion, which could in principle be due 
to scale-breaking terms in the pressure or residual external potentials, this is a 
genuine dissipative effect, since vorticity is conserved in ideal hydrodynamics. 

The main difficulty in extracting the shear viscosity from the analysis of scaling 
flows is associated with the role of the dilute corona of the cloud. Kinetic theory 
predicts that in the dilute limit the shear viscosity is independent of density and only 
depends on temperature. A simple analysis of the type presented in Sec. 13.31 then 
implies that the dilute corona does not generate a dissipative force. It nevertheless 
dissipates a large amount of energy. The analysis also suggests that freezeout only 
occurs very late, see Sec. 13.41 There are a number of aspects of this analysis that 
need to be improved: 
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• The Navier-Stokes equation is based on the assumption that the viscous correc- 
tion to the stress tensor appears instantaneously. This is particularly problematic 
in the case of scaling flows, because the viscous contribution is spatially constant. 
The fact that the ideal stresses propagate outward with the expansion of the sys- 
tem whereas the dissipative stresses appear immediately indicates that causality 
is violated. This problem can be addressed by including a finite relaxation time, 
or by solving a more complete set of second order hydrodynamic equations. 

• We have studied the effect of dissipative forces in the Navier-Stokes equations, 
but we have computed the non-dissipative forces (pressure gradients) based on 
an approximately isentropic expansion. This procedure neglects reheating, and 
violates energy conservation. Reheating is important in the dilute corona, and 
breaks the scaling nature of the expansion. 

In addition to implementing these technical improvements it is important to con- 
sider other experimental setups that are directly sensitive to the spatially constant 
part of the shear viscosity. One option would be to measure the attenuation of sound 
propagating in a very long elongated trap. Another idea would be to directly mea- 
sure the decay of a shear flow in a long channel. 

Finally, we summarize the existing experimental constraints on the shear viscos- 
ity of the unitary Fermi gas: 

• The damping of collective oscillations constrains the trap average (7])/5 = (ct s ). 
We find that this quantity varies between (a s ) ~ 1 at T/Tp ~ 0.8 and (a s ) ~ 0.5 
at T/Tp ~ 0.2. In the regime 0.4 < T/Tp < 0.8 the temperature dependence is 
consistent with rj ~ (ml) 3 / 2 and a relaxation time that scales as Tr ~ T]/(nT). 
At lower temperatures an additional contribution is needed. In a simple model 
the minimum of the shear viscosity to density ratio is rj/n ~ 0.2. 

• The expansion of a rotating cloud gives (a s ) ~ 0.8 at T/Tp ~ 0.8, and (a s ) ~ 
(0.0 — 0.4) at T/Tp ~ 0.2 1321 . The latter results are smaller than the values ex- 
tracted from collective oscillations, although the errors are also somewhat larger. 
It will be important to determine whether this discrepancy is due to the effects of 
the dilute corona, and whether the smaller values of (a s ) are more representative 
of the shear viscosity to entropy density ratio in the core. 

Note added: After the initial version of this contribution was finished dissipative 
effects in the expansion of a dilute Fermi gas at temperatures T 3> Tp were stud- 
ied experimentally by Cao et al. [l40l . This work nicely demonstrates the scaling 
(o£„) ~ T 3 predicted in Fig. [4] Numerical solutions to the equations of dissipative 
hydrodynamics were studied in BT1 . This work shows that quantitative estimates of 
the shear viscosity have to take into account the effects of reheating. 

Acknowledgments: This work was supported in parts by the US Department of 
Energy grant DE-FG02-03ER41260. We are grateful to John Thomas for many use- 
ful discussions, and to Jiunn-Wei Chen for pointing out an error in an earlier version 
of this contribution. 



Scaling Flows and Dissipation in the Dilute Fermi Gas at Unitarity 



31 



References 



1. I. Bloch, J. Dalibard, W Zwerger Rev. Mod. Phys. 80, 885 (2008) |arXiv:0704.301 1 1. 

2. S. Giorgini, L. P. Pitaevskii, S. Stringari, Rev. Mod. Phys. 80 12 15 (2008) |j rXiv:0706.3360|. 

3. T. Schafer and D. Teaney, Rept. Prog. Phys. 72, 126001 (2009) I arXiv:0904 3l07l rhep-phll. 

4. P. Danielewicz and M. Gyulassy, Phys. Rev. D 31, 53 (1985). 

5. G. Policastro, D. T. Son and A. O. Starinets, Phys. Rev. Lett. 87, 081601 (2001) 
|arXiv;hep-th/010406"6) . 

6. P. Kovtun, D. T. So n and A. O. Starinets, Phys. Rev. Lett. 94, 111601 (2005) 
(arXiv:hep-th/040523 f) . 

7. K. M. O'Hara, S. L. Hemm er, M. E. Gehm, S. R. Granade, J. E. Thomas, Science 298, 2179 
(2002) |cond-mat/0212463l . 

8. J. Kinast, S. L. Hemmer, M. E. Gehm, A. Turlapov, and J. E. Thomas, Phys. Rev. Lett. 92, 
150402 (2004) lcond-mat/04035401 . 

9. J. Kinast, A. Turlapov, J. E. Thomas, Phys. Rev. A 70, 051401(R) (2004) lcond-ma t/0408634 1 . 

10. J. Kinast, A. Turlapov, J. E. Thomas, Phys. Rev. Lett. 94, 170404 (2005) |cond-ma t/0502507 |. 

11. A. Altmeyer, S. Riedl, C. Kohstall, M. Wright, R. Geursen, M. Bartenstein, C. Chin, J. Hecker 
Denschlag, R. Grimm, Phys. Rev. Lett. 98 040401 (2007) lcond-mat/06093901 . 

12. A. Altmeyer, S. Riedl, C. Kohstall, M. Wright, J. Hecker Denschlag, R. Grimm, Phys. Rev. 
Lett. 98 103602 (2007) lcond-mat/06 112851 . 

13. M. J. Wright, S. Riedl, A. Altmeyer, C. Kohstall, E. R. Sanchez Guajardo, J. Hecker Den- 
schlag, R. Grimm, Phys. Rev. Lett. 99, 150403 (2007) |arXiv:0707.3593rcond- mat.other]]. 

14. B. Clancy, L. Luo, J. E. Thomas Phys. Rev. Lett. 99 140401 (2007) larXiv;0705.2782i llcond- 
mat.other]]. 

15. S. Riedl, E. R. Sanchez Guajardo, C. Kohstall, A. Altmeyer, M. J. Wright, J. Hecker 
Denschlag, R. G rimm, G. M. Bruun, H. Smith, Phys. Rev. A 78, 053609 (2008) 
|arXiv:0809. 1 8 14 r cond-mat.otherl 1 . 

16. J. Carlson, S. Reddy, Phys. Rev. Lett. 95 060401 (2005) lcond-mat/05032561 . 

17. A. Bulgac, J. E. Drut and P. Magierski, Phys. Rev. A 78, 023625 (2008) larXiv:0803.3238l 
[cond-mat.stat-mech]]. 

18. S. Nascimbene, N. Navon, K. Jiang, F. Chevy, C Salomon, |arXiv: 091 1.0747 cond-mat.quant- 
gas]]. 

19. C. Menotti, P. Pedri, S. Stringari, Phys. Rev. Lett. 89, 250402 (2002) |cond-mat/0208150 1. 

20. J. E. Thomas, J. Kinast, A. Turlapov, Phys. Re v. Lett. 95, 120402 (2005) lcond-ma t/0503620 1. 

21. D. T. Son, Phys. Rev. Lett. 98, 020604 (2007) larXiv:cond-mat/05 1 1 72 1 1 . 

22. L. Luo, J. E. Thomas, J. Low Temp. Phys. 154, 1 (2009), jiarXiv:081 1.1 159 cond-mat. other]]. 

23. G. Rupak and T. Schafer, Phys. Rev. A 76, 053607 (2007) larXiv;0707. 15201 [cond- 
mat.other]]. 

24. G. M. Bruun, H. Smith, Phys. Rev. A 72, 043605 (2005) lcond-mat/05047341 . 

25. G. M. Bruun, H. Smith, Phys. Rev. A 75, 043612 (2007) |cond-mat/0612460| . 

26. L. S. Garcia-Colina, R. M. Velascoa, F. J. Uribea, Phys . Rep. 465 149 (200 8). 

27. T. Nikuni, A. Griffin, Phys. Rev. A 69, 023604 (2004) |cond-mat/ 0309269 |. 

28. A. Griffin, T. Nikuni, E. Zaremba, "Bose-condensed gases at finite temperature", Cambridge 
University Press (2009). 

29. P. Romatschke, ||arXiv;0902.3663i rhep-phl. 

30. G. M. Bruun, H. Smith, Phys. Rev. A 76, 045602 (2007) |arXiv:0709. 16171 . 

31. M. Edwards, C. W. Clark, P. Pedri, L. Pitaevskii, S. Stringari, Phys. Rev. Lett. 88 070405 
(2002). 

32. J. E. Thomas, |arXiv;0907.0140V 2 [cond-mat.quant-gas]. 

33. B. Clancy, Ph.D. Thesis, Duke University (2008). 

34. H. Heiselberg, Phys. Rev. Lett. 93, 040402 ( 2004) lcond-mat/040304fl : 

35. S. Stringari, Europhys. Lett. 65, 749 (2004) I cond-mat/03 1 26141 . 

36. A. Bulgac and G. F. Bertsch, Phys. Rev. Let t. 94, 070401 (2005) lcon d-mat/04046871 . 

37. T. Schafer, Phys. Rev. A 76, 063618 (2007) larXiv:cond-mat/0701 25 1 1 . 



32 



T. Schafer and C. Chafin 



38. B . A. Gelman, E. V. Shuryak, and I. Zahed, Phys. Rev. A 72, 043601 (2005) 
lnucl-th/04100671 . 

39. A. Turlapov, J. Kinast, B. Clancy, L. Luo, J. Joseph, J. E. Thomas, J. Low Temp. Phys. 150, 
567 (2008) larXiv:0707.2574l . 

40. C. Cao, E. Elliott, J. Joseph, H. Wu, J. Petricka, T. Schafer and J. E. Thomas. larXiv: 1007.26251 
fcond-mat.quant-gas] . 

41. T. Schafer, arXiv: 1008.3876 [cond-mat.quant-gas]. 



